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O ,' Abstract 

D ■ 

f^ , ' We discuss a specific cut-ofi' efi'ect which appears in applying the non- 

perturbative RI/MOM scheme to compute the renormalization constants. 
To illustrate the problem a Dirac operator satisfying the Ginsparg- Wilson 
^ ' relation is used, but the arguments are more general. We propose a simple 

("V^ , modification of the method which gets rid of the corresponding discretiza- 

C^ ■ tion error. Applying this to full-QCD simulations done at a = 0.13 fm 

with the Fixed Point action we find that the renormalization constants 
are strongly distorted by the artefacts discussed. We consider also the 
f~^ ' role of global gauge transformations, a freedom which still remains after 

c7^ ' the conventional gauge fixing procedure is applied. 
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1 Introduction 



The correlators of the bare fields and composite operators in lattice simula- 
tions are not completely physical. Their amplitudes depend on the details of 
the action used in the simulation and from the actual form of the operators 
involved. To be able to compare different simulations (with different lattice 
spacings, actions, operators, etc.) and also to compare the numerical results 
with perturbation theory one needs some convention to bring these quantities 
to a "common denominator" . For that purpose one has to introduce renormal- 
ization constants, and define the renormalized fields and operators. These will 
not depend any more on the microscopic details (apart from the discretization 
errors) but will depend on how the renormalization conditions have been chosen, 
i.e. on the renormalization scheme. 

In continuum perturbation theory one generally uses dimensional regular- 
ization and the MS scheme to fix the convention of removing (factoring out) 
the divergencies. In lattice simulations this method is not applicable, and to be 
able to connect to the MS scheme one has to use some condition to define the 
renormalized fields in a way which is independent of the actual regularization. 
Such scheme is the non-perturbative Rome-Southampton RI/MOM technique 
[1]. In this scheme one fixes some renormalized Green's functions to their free 
theory counterpart at some scale p^ = /^^ where p is the momentum of external 
particles and /i is an appropriately chosen scale. QCD has an extra complication 
due to gauge invariance: the quark correlators are not gauge invariant, hence 
they average to zero when integrated over different gauge-equivalent copies. 
Therefore in the RI/MOM scheme one uses a gauge fixing condition to make 
sense of the Green's functions with external quark legs. This follows closely the 
perturbation theory where one has to fix the gauge anyhow. 

The RI/MOM scheme has become the standard way of fixing the renor- 
malization constants, and has an extended literature (see [21 El [H IH [6l [7] and 
references therein). In particular, there are many versions to formulate the same 
renormalization condition. These coincide in the continuum limit but contain 
different lattice artefacts. 

We believe we found a specific artefact which can have a large effect on the 
renormalization constant and (as far as we know) all actual definitions used in 
the literature suffer from it. Fortunately, one can cure the problem in a simple 
way. We discuss the issue for the case of the overlap Dirac operator where this 
artefact is most transparent. However, the problem (and the cure) is general. 

The outline of the paper is the following: After stating the conventions we 
describe the Landau gauge fixing, and suggest a simplification related to the 
remaining global gauge degrees of freedom. 



The artefact in the renormahzation constants is discussed in detail on the 
example of the quark field and the composite bilinear operators. 

Having modified the matching relations, we compare different choices for the 
renormahzation constants of the quark field and of the covariant scalar density 
in an actual simulation done on a 12* lattice with 2+1 light flavors with the 
Fixed Point (FP) action. 



2 Renormalization 



We use the following conventions for the renormalization constants. The renor- 
malized fermion field is related to the bare field aqj 

while the renormalization factor for local operators enters aqj 

O^^ZoO. (2) 

To implement the RI/MOM scheme a gauge fixing is necessary. One averages 
then the propagators over the gauge-fixed configurationqj. The fermion propa- 
gator on a given gauge field configuration is denoted by S{U) and its gauge-field 
average by 

S - {SiU)} (3) 

According to eq. ^ this is renormalized as 

S^ = Z-^S. (4) 

2.1 Gauge Fixing and color averaging 

Eq. ^ and other gauge-field averaged non-gauge-invariant correlation functions 
considered below assume gauge fixing. We use the Landau gauge which is de- 
fined by maximizing the sum 

tTj2K{^) + Ul{x)) (5) 

with respect to gauge transformations. Here (/^{x) = e^^^'^^^ £ SU(3) is the 
lattice gauge field variable. In the continuum limit this condition is reduced to 
d^A^ix) = 0. 



^ This convention agrees with that of different textbooks [8] [Q] 1101 1111 

^This convention is widely followed in the general and lattice literature (an exception is [9]) 

•^generated according to the appropriate Boltzmann weight 



A global gauge transformation C/^(x) -^ ^U^j,{x) = gU^{x)g\ where g G 
SU(3) is independent of x, rotates the gauge fields while keeping the value of 
the functional in eq. © unchanged. The extent of the effective averaging in this 
global degree of freedom by the available configurations depends not only on the 
number of configurations, but also on the actual algorithm for fixing the gauge. 
It is useful and natural to integrate over g explicitly. Consider, as an example, 
the propagator S{x, y; U) on a given configuration. Performing a global gauge 
transformation we get 

S{x, y- U) -. S{x, y; ^U) = gS{x, y; U)g^ . (6) 

When we are interested in S'(x, y), which is the average of the individual propa- 
gators over the gauge fields in Landau gauge, as argued above, we can addition- 
ally average this quantity over the global gauge freedom. This "color averaging" 
amounts to replacing the color sub-matrices (at given Dirac indices) by 3 x 3 
matrices which are trivial in color, according to 

S{x, y) -^ dg gS{x, y)g'^ ^ ^^ ^ *''= '^(^' ^^ ' ^'^^ 

where Ic is the unit color matrix and trc denotes trace in color indices. 

In the following we shall perform an additional color averaging whenever an 
averaging over gauge fields is done. In other words the gauge-field averaged 
Green's functions will be considered to be 4 x 4 matrices having only Dirac 
indices. 

Note that to perform a color averaging is mainly a matter of convenience, it 
is easier to work with 4x4 than with 12 x 12 matrices. In practice the actual 
color sub- matrices are already nearly proportional to Ic. 



2.2 The renormalization conditions 

The RI/MOM technique [1] consists of imposing the condition that some renor- 
malized Green's functions at a given scale p^ — fi^ are equal to their correspond- 
ing tree level values. 

One would like to connect this renormalization scheme to other schemes by 
using perturbation theory and, at the same time, control the cut-off effects on 
the lattice. These conditions constrain the scale ^: 

M -^ H-^-, (8) 

a 

where M is a typical non-Goldstonc boson mass scale in QCD. 



We can formulate an analogous condition in configuration space by request- 
ing that the renormalized Green's functions should be matched at distances 
|a; — 2/|, where 

a< [x- j/| = - < Ifm. (9) 

A* 

Since in full QCD applications eq. ^ is satisfied only marginally one needs 
further decisions in the applications in order to avoid large cut-off effects. 



2.3 A special cut-ofF effect rooted in lattice chiral symme- 
try 

2.3.1 The quark field renormalization constant Z' in the RI' scheme. 



We shall illustrate this cut-off effect and the way to eliminate it in the context 
of the quark field renormalization factor Z' which is defined in the so called RI 
schemeQ In this case the Green's function considered is just the propagator in 
Fourier space given by 

x,y x,y 

The generalization to Green's functions with a bilinear operator will be treated 
the next section. 

Following the general procedure of matching with the free theory one might 
request 

5^b)|p.=^. = z'-' s{p\.^^. .. s{py^^%^^, , (11) 

where S^'^"'^ is the inverse of the Dirac operator on the trivial gauge configuration 
U = I and /i^ is constrained by eq. ([8]). One might also request the same 
equation in coordinate space using eq. ([9]). Eq. (fTTjl is a relation between two 
4x4 matrices. If n^ satisfies eq. ([8]) and the statistical error in the simulation 
is small, one might expect that all elements of this matrix equation can be 
(approximately) matched by a single parameter Z' . This, however, does not 
hold in general. 

The symmetries of the lattice action and the gauge fixing condition imply 
that only the 1 and 7^ Clifford algebra elements enter in the propagator 

S{p)^bo{p)l + ib^{ph^. (12) 



^In the RI' scheme one uses S{p) ^ while the RI scheme uses dS{p) ^ /dpf^ in the matching 
condition. 



The presence of the unit 4x4 matrix 1 is related to the t/(l) anomaly and the 
problem of doublers [12] . 

Consider a Dirac operator which satisfies the Ginsparg- Wilson (GW) rela- 
tion [H M 

l5S{U)^y + S{U)^y-fz = 2ai?(L/),y75 , (13) 

where R(U)xy is a local operator which is trivial in Dirac indices. The lattice 
spacing a is written out here explicitly to indicate that this term vanishes in 
the formal continuum limit. In particular, consider the simplest case given by 
Neuberger's overlap Dirac operator [TS], 

R{U)^y = Kd^y . (14) 

where k is a real number of order 1. These equations imply that the coefficient 
of the unit Dirac matrix 1 in the propagator D~^(U)xy — S(U)xy is anSxy 
Averaging over the gauge configurations U in Landau gauge gives 

IbSxy + Sxy^5 = 2aKSxy^5 ■ (15) 

The expectation value S^y depends only on x — j/ in a periodic box. We get for 
the coefficient of 1 in coordinate and momentum space, 

boix^y) ^ anSxy (16) 

and 

bo{p)^aK, (17) 

respectively. The matrix S{py'^'^'^ has the same Dirac matrix structure as S{p) 

sipf = h^ipt^n + ih^ipt^^^^ , (18) 

with 

boipY'^'' = boip) = aK , Vp. (19) 

Matching the 6o part in eq. pT|) would give Z' = 1, which is obviously nonsense. 
Eqs. (dni) and P?|) give the explanation: bo is 0(a), a pure cut-off effect in this 
context. In addition, d ^ \x — y\ = violates eq. ^. Obviously, the 1 part of 
the propagator, whose presence is absolutely essential to avoid doublers and get 
exact chiral symmetry, should not be taken (alone, or in combinations with the 
7^ part) in determining Z' . For this reason the conventional matching condition 
proposed in jTj and widely used in the literature, 



z'^^tv{s{p)-Hj,b,{py^^^) 



= 1 (20) 



^2 ^((2 



is plagued by this artefact since the 7^ part of S{p) ^, unlike the 7^ part of 
S{p), contains bo{p). The same applies to an alternative definition 

Z;tr(5(p)-i5(p)^'--)| .^. = 1 (21) 



proposed in [3] 

Although we treat cut-off effects here, these distortions, as we illustrate 
below, can be large if the lattice unit a is not very small. 

In the arguments above we discussed the case of the overlap operator. In this 
case the considerations are fully transparent. The message is, however, general: 
the 1 part of the propagator S{p) is a special cut-off effect which should not 
enter the renormalization conditions. 

To illustrate the point we use the results from the simulations with 2-f 1 light 
fermions on a 12'' lattice, using the Fixed Point (FP) Dirac operator at lattice 
spacing a = 0.13 fm [16l[T7] 



We plot in Fig.mthe ratios 6o(p)/6ob)^'"' and b ^{p) / b ^(j))f'''"' for /i = 1, . . . , 4 
at different momenta p — ipi,P2,P3,P4) (not averaged over the directions of 
p). While the ratios for /i — 1, ... ,4 are consistent with each other, the &n 
ratio differs significantly from the others, in agreement with our expectationsO 
The 4-component vectors bfj,{p) and &^ (p)-'^^'^'^ for a given p are parallel to each 
other to a surprisingly high accuracy: for the angle 9 between them we find 
1 — cos9{p) < 10~^ which corresponds to an angle 9 < 0.3°. 

Due to this fact the matching conditions for the 7^ parts are consistent with 
each other. To be specific, we propose to average the four equations bfi{p)/Z'^ = 
b^u'^'^ip) with the weights b^{pY^^'^ , i.e. to use the renormalization condition 






Eb^ripWrip) 

Here 



(22) 



b^{p) = ^tr{S{ph^). (23) 

It is convenient to introduce the operation of subtracting the unit Dirac 
matrix pariij 

M^M- Q tro m\ 1 , (24) 

With this notation S{p) — ib^{p)^^ and eq. ([^^ can be written as 

Z',= t^{^{p)S{p)j;^)l^^^^. (25) 

^While our FP Dirac operator satisfies the GW equation to a good accuracy, the 1 part 
of the propagator is a nontrivial (nearly local) operator R{x, y; U) hence feo(p) is expected to 
depend on p. 

^For later application we allow the matrix M to have color indices as well. The color part 
is not affected by this procedure. 




p [GeV": 



Figure 1: Different coefficients of the propagator S{p) — bo{p)l + ibu{p)^v 
divided by tfie corresponding coefficients of the free propagator. 



For consistency it is useful to check the other choice (obtained by using h^{p) 
as weights) 

l-=ir{S{p)-^^{p)free\,^^^. ■ (26) 

The two definitions should coincide if hf^{p) is indeed parallel to b^{p)free- 

Fig. [5] shows the results for Z' obtained by different definitions, eqs. ([^5]) . 
(EOD and dm). 



2.3.2 Eliminating the special cut-off effects in ZrZq 



We shall consider bilinear operators of the typqj 

Or(x)=^^yr(a:;[/)j,,V., 



(27) 



where the kernel r(a;; U) is local, its value is essentially zero if the distances 
I a; — 2/1 or |x — z| are larger than 0(a). (For the naive densities and currents it 
has the form r(x; U)yz = TSxySxz, where F is the corresponding Clifford algebra 
element.) 



^We assume for simplicity that O is non-singlct in flavor. 
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Figure 2: The values of Z' calculated by different methods: from cq. ((25)) (thick 
line), and by two conventional methods, given in eq. ([20]) (thin line) and eq. (j2ip 
(dashed line). 



The Green's function we consider i^l 



Gv{v) = ^E^' 



u,v 

y,z 



i(u—v)p 

V 



{ipu^v^{x;U)y^il!z')py) 



= ^E^"^^""'^'^ T.(^iU)uyrix; U)y.SiUU) (28) 

When one considers matching with the corresponding free field value in coordi- 
nate space it is obvious that the distances \u — y\ and |z — w| should be taken to 
be physical, much larger than a. As discussed previously, performing the match- 
ing in momentum space also contains the non-physical contact term okI in the 
case of the overlap and some similar nonphysical contribution for other Dirac 
operators. It is natural to omit these contributions in the matching condition. 



Using the notation (j24p for subtracting the 1 part of the propagator we 
define the modified Green's function 

^r(p) - ^T.^-'^"'^^'^ Y.(^iUUT{x; [/),.^(t/)..) . (29) 



V' 



y,z 



V 



*Note that S{U)xy is still a 12 X 12 matrix; the color averaging is done together with 
averaging the whole product over the gauge configurations. 



It is convenient to define the corresponding amputated Green's function 

Ar{p) ^^{pr'Gr{p)S{pr\ (30) 

wliere we use S{p) for amputation by the same argument as discussed above. 
This is the quantity which we shall use in the renormalization condition. We 
require that the matrix equation 

ZrZ.Arip) ^ ^ ^ Mpf'' , , (31) 

holds approximately. 

Correspondingly, one can use a scalar renormalization condition 

ZrZ,tr(rAr{p)) ^ = tr (rAr(p)^'--) ^ ^ (32) 



for the product Z^Zq. In other words, we replace Ar(p) in the conventional 
lattice definition of Z^Zg by Ar(p). 

To illustrate the cutoff effect discussed here consider the covariant scalar 
operator for a general GW Dirac operator. (Our parametrized FP operator is 
an approximation to such Dirac operator.) In this case one has 

S^yiU) = R,yiU) +S^y{U) , (33) 

where R is a, local operator, proportional to 1 and appearing on the r.h.s. of the 
GW relation eq. p3)) . while Sxy{U) is traceless in the Dirac indices. The covari- 
ant scalar operator is given essentially ITS' by Os = ^'I'V' with T = 1/{2R) k, 1. 
Inserting this into eq. (j28p one obtains four terms, symbolically written as 
Gs = (StS) + {RTR) + cross terms = G + G^^^ + . . .. For the proper matching 
one needs only the first term, the second should be omitted. (The cross terms 
have "wrong" Dirac structure and do not show up in the actual matching.) 
Roughly speaking, Gs ~ -R/2 — Ip /{2R) « 1/4 - b"^, where S ~ ih^^^. Here one 
expects that R(p) is approximately constant while h'^{p) ^ 1/p^ decreases with 
increasing p. 

In Fig. [3] we plot the full Gs and the terms Gs and Gg using our data. The 
figure shows that the nonphysical part Gg dominates for p^ > 5 GeV'^ hence 
taking the full quantity Gs instead of Gs results in a very strong distortion of 
the true matching ratio. Fig. 2] illustrates how the proposed modifications affect 
the value Zi-Zq. 
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Figure 3: The full unamputated Green's function G{p) ^ (STS) (thin line), its 
unphysical part G''-"'^-'(p) ~ (RTR) (dashed line), and the physical part Gr{p) ~ 
(STS), for the scalar operator, F = 1. (thick line). 



2.3.3 The quark field renormalization constant Zq in the RI scheme. 



Consider now the special case when the bilinear operator Or{x) is the conserved 
current (related to the flavor conservation): 



Vnix) =tpyK^{x;U)y^ip^ 



(34) 



Given the Dirac operator D{U), there is a simple procedure to find Kfj^ [TO1IT5]. 
The kernel satisfies the Ward identity in coordinate space 



Y^ d;K^{x- U)y. = {5,y - 4.) D{U) 



yz , 



(35) 



where d* is the backward lattice derivative in x. Starting from eq. ([55)1 it is 
not difficult to get the Ward identity for the corresponding Green's function in 
momentum space 

^''^'^P^^~'W,^^^'^' (36) 

where G^''{p) is defined in eq. ([28]) with T{x)yz -^ K^{x; U)yz. For the ampu- 
tated Green's function we get 



d 



K^ip) - S{p)-'G}:^{p)Sip)-' = ^^S{p) 



(37) 
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Figure 4: The values of ZsZq calculated from eq. 
conventional method (thin line) . 



(thick line), and by the 



The quantity on the r.h.s., dS{p)^^ /dpf^ is used to fix the rcnormalization fac- 
tor of the quark field. The corresponding condition defines a scheme which is 
different from that given by eqs. (fTTj) and ((22| . The quark field rcnormalization 
factor in this scheme is denoted by Zq. The Ward identity shows that in this 
scheme the conserved current does not renormalize, i.e. Zy^ = 1- (The advan- 
tage of using the vertex function of the conserved current lies in the fact that it 
avoids approximating the derivative over the momentum by discrete derivatives 
available on a finite lattice.) 



Note that the conventional definition of Zq through eq. ([37] 



Zq tr 7^ 



dSipY 



dptj 



= tr 7m- 



dS{p)7ri 



free 



p^—f_L^ 



dpf. 



(38) 



also suffers from the same lattice artefact. Instead of this condition one should 
take _ 

^ dS{p)-' dS{j>)ll 



' dp. 



dpn 



P'^—f_L'^ 



or equivalently 

Zq^ipr^G]:{p)sipy 



S{p)free ^ ^.^ {p) free S {j)) f^l^ 



(39) 



(40) 
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To have a scalar equation it is convenient to take out the 7^ part and sum over 



Zq tr 



(7^5(p)-^<^(p)5(p)-i 



= tr 



\lt^S[p)^,^^fi^; ip)freeS{p)f^l 



(41) 



T=7Vc 



Note that here we have G^" which at first sight differs fi-om Gjj^ defined 



in eq. (I29p . However, due to the Ward identity the Green's function for the 
conserved current G^{p) is effectively linear in the propagator S{p) and the 
two quantities coincide. For general operators, however, one should use the 
quantities Gr or Ap. 




Figure 5: The values of Zq calculated by different methods: from eq. (j40| (thick 
line), and by the conventional method (thin line). 



3 Summary 



We suggested to color average the gauge-field averaged propagators and Green's 
functions used in the RI/MOM scheme, reducing them to 4 x 4 matrices having 
only Dirac indices. This is, however, just a matter of convenience. 

Further, we argued that the presence of the unphysical, short range part of 
the quark propagator proportional to the unit 4x4 matrix (in Dirac indices) 
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is responsible for a considerable 0(a) artefact in the renormalization constants 
of the quark field and the composite operators. We illustrated that by omitting 
this contribution one can satisfy the matching condition in the matrix sense to 
a remarkable precision, which is spoiled otherwise by the non-physical contri- 
bution. 

The effect discussed here is quite large for the lattice simulations with the 
Fixed Point fermions at a lattice spacing a = 0.13 fm. However, the discretiza- 
tion error is expected to be less pronounced in simulations at smaller lattice 
spacing. It can also be smaller for the Wilson action where the coefficient of the 
1 in the propagator can decay faster in p. These questions need further study. 

Acknowledgement We thank Peter Hasenfratz for enlightening discussions 
on the problem of these artefacts, Anna Hasenfratz for sharing her experience 
with the Rl/MOM scheme and Dieter Hicrl for helping us with the numerical 
data. 
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